Principles of statistical mechanics: the energy duality 
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We argue that statistical mechanics of systems with relaxation implies breaking the energy func- 
tion of systems into two having different transformation rules. With this duality the energy approach 
incorporates the generalized vortex forces. We show general theorems and their implications and 
apply the approach to the particle confinement in fields of rotational symmetry. Misconceptions of 
extensive use of the quasienergy and generalized thermodynamic potential theories are exposed. 
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Introduction 

While the main word in physics is interaction of sys- 
tems given by an energy function, this notion is well de- 
fined for the ideal of closed systems via their motion inte- 
grals. For the real world, of open systems in long lasting 
non-equilibrium states against the background of fluc- 
tuating environment, the meaning of energy function is 
vague and causes misleading associations and mistakes 
in energy reasoning and modeling. A new formulation 
of energy approach to open systems is suggested in this 
Letter. 

We found this matter important even in conditions 
of low-rate irreversible forcing, for its nonlinear cumu- 
lative effect on the trend of system's steady state, sta- 
bility and fluctuations can be drastic and different from 
the trend given by the theory of generalized thermody- 
namic potential [1-3] commonly accepted in the study of 
heating/cooling, transport and phase transitions. This 
paradigm fits the ideal of closed systems, showing no en- 
ergy vagueness. The generalized potential of a system 
relaxing in steady conditions to a probability density p s t 
connects to p s t by 



Pst 



Ne 



N~ 



>dY 



(1) 



where the integral is over the volume F of system phase 
space and the reversible motion is on surfaces $ = const. 
The properties of the system mainly depend then on the 
local properties of the minima of The analogous ap- 
proach to systems under high frequency fields is in terms 
of the picture where the hf field looks fixed or its effect 
is time-averaged. In all this, Eq. (1) merely redefines the 
distribution p st in terms of function <i>, and taking this 
function as the energy integral of reversible motion pro- 
vides the physical basis of the theory, but implies rigid 
constraints. 

Let us start the analysis of these constraints and the 
energy approach beyond them with formulating general 
theorems. 
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Three theorems: of entrainment, energy 
indeterminism and energy duality 

I. The system whose reversible- motion integral and 
steady distribution of states p s t are in a one-to-one re- 
lation must carry along, on the average over p st , the 
environment contributing to its fluctuations and dissi- 
pation. The relation of p st to $ can be taken as the one 
in question without loss of generality. Connecting $ to 
the energy function of system implies scaling this func- 
tion in terms of environmental-noise level of energy. The 
energy scale set so must vary proportionally with the en- 
ergy in arbitrarily moving frames while $ defined by (1) 
must not vary. This constraint can hold only for the sys- 
tems carried along with the environment on the average 
and must break beyond this entrainment ideal, hence the 
proof. 

II. The energy integral of reversible motion of the sys- 
tems relaxing to a steady distribution p st beyond the en- 
trainment ideal ceases to exist. Indeed, as follows from 
the proof of theorem I, neither $ (1) nor any univalent 
function of it can be the energy integral, and it cannot 
be a function of transient process to p s t as well, for the 
energy integral of reversible motion must be independent 
of transient states in Y. 

The energy integral break-up can be proved also as 
follows. The non-entrained steady state implies a steady 
mean motion in Y. The relaxation to it implies that the 
irreversible forces exerted on the motion do not vanish as 
t — ► oo and are of vortex type in Y as their forcing toward 
the steady motion and against it differ in sign. Such 
forces form both the p st and the approach to it. Thus 
the energy integral ceases to exist there, being blurred 
by the vortex forces. 

III. The whole blurring must come down to that the 
notion of system energy function (Hamiltonian) adopted 
as one- valued for a certain standing must then bifurcate. 
Its two Hamiltonians represent the constraints relating 
the partial derivatives of Hamiltonian with respect to the 
phase space variables of the system to its instant state. 
One is bound with <E> and represents the constraints on 
the motion relative to the environment as the source of 
diffusion/dissipation; the time reversal decomposition is 
with respect to the parity associated with the partial 
derivatives of just this Hamiltonian, so the irreversible 
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drift is determined by them and hence this Hamiltonian. 
The other governs the system motion unrelated to the 
said source, but is also involved in forming the p s t and 
the transition to it. This is the only way to comply with 
the arguments in I and II, and this is what we called the 
energy duality. 

The difference between the two energy functions is due 
to the work of vortex forces. In its turn the difference 
function acts as their energy measure. Thereby the vor- 
tex forces are incorporated into the energy concept. 

The theorems of entrainment, indeterminism and du- 
ality shown above admit extension to the systems in un- 
steady conditions so long as the probability description 
is adequate and the probability distribution of system 
states exists and relaxes from the domain of initial con- 
ditions to a common, limit i-dependent distribution. The 
proof follows from the fact that a limit distribution is re- 
ducible to a steady one with univalent transformations 

of r. 

Historically (e.g. reviews [3,4]) the necessary and suffi- 
cient conditions of generalized potential theory were for- 
mulated as detailed balance within the framework of au- 
tonomous Fokker-Planck equations for the variables and 
parameters dividable into odd and even with respect to 
time reversal. The detailed balance complies with the en- 
trainment ideal. Our theorems presented above set the 
energy approach to wider conditions. 

Let us now give the introduced notions substance and 
consider implications. 



Problem statement via kinetic equations 

We consider the probability density p(z,t) of system 
states z as prime and determine the concept of system 
energy function from its time evolution. The interactions 
with fluctuating environment will be considered within 
the equations of general form 



dp , , dSi 



(2) 



from a given initial distribution p(z, 0) under natural 
boundary conditions. Without the part dSi/dzi, Eq. (2) 
is the Liouville's equation of a dynamic system of Hamil- 
tonian H , a function of a number of canonical coordinates 
x = (a?i, . . . , x n ) and moments p — (pi, . . . ,p n ) taken as 
z, z = (x,p). [, ] denotes the Poisson brackets, 



[H,p} = 



3H_dp_ _ dH dp 

dxi dpi dpi dxi ' 



Summation over repeated dummy indexes is implied. 
S = {Si} is a 2n-vector functional of p in z, vanish- 
ing at the boundaries to preserve the normalization of 
p, includes the irreversible probability currents. These 
currents depend on the system's motion, and this depen- 
dence in its turn affects the energy concept. 



In diffusional approximation 
d 



S i =[f i -d ik — U , 



i 7 k = 1, . . . 2n 



(3) 



where / = {/,} are dissipativc forces and d = {dik} is 
a symmetric positive semi-definite matrix of diffusion. 
Eq. (2) is then a general Fokker-Planck equation written 
in terms of canonical variables of the reversible dynamics 
given by H. 

S may be nonlinear in p and integro-differential in z, 
unless specified otherwise, and we use the term kinetic 
also for such equations (2). The parameters entering H 
and S may depend not only on z but also on t to allow 
for interactions with varying regular fields and varying 
chaotic environments. 

Entrainment ideal. Its simplest modeling corresponds 
to S of form (3) with / and d having non-vanishing com- 
ponents pertaining only to the space of p, with d a con- 
stant positive-definite matrix n x n and / a n-component 
force of viscous friction 



fi = -PdikVk 



(4) 



where (3 1 is a scale of noise energy and v = (v\, . . . ,v n ) 
the velocity given by Hamiltonian dynamics, 



Vi 



dH(x,p,t)/dpi. 



(5) 



For a ^-independent H bounded below, the stationary 
solution of (2) is then the Boltzmann distribution 



p B (z) = Ne~^ z \ N~ x =Ji 



-PH(z) d 2 



(6) 



The meaning of f3 expounds the known equipartition the- 
orem ensuing from (6): for every component of z (coor- 
dinate or momentum) whose contribution to H reduces 
to a square term, say, Ki(zj — K2) 2 with k\ 7^ and 
independent of Zj, its mean over the Boltzmann distri- 
bution comes to 



(Ki(Zj - k 2 ) 2 } b =0- 



(7) 



irrespective of parameters n\^- This theorem holds not 
only for constant but also for n\ t 2 depending on other 
components of z and on t, which is not insignificant for 
"quasistationary" Boltzmann distributions. By virtue of 
(5) and assuming natural boundary conditions 



J ve -/3H( x ,p,t) d n p 



0. 



(8) 



So the average (v(x,p, t)) B = over p for any x, which is 
the entrainment ideal for the case. 

The same features hold beyond the linear friction and 
constant diffusion rates, being also in effect for the / 
and d of elements pertaining to both x and p spaces and 
depending on both z and t, provided that 



fi = -f3du 



dH 

dzk 



i,k = 1, 



,2n 



(9) 
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and that the matrix d ensures approach of the solutions of 

(2) to a unique steady state. One comes then to (6) and 
thus to the proof. Note that the 2n vector v = dH/dz 
has norm \v\ = \[z,H]\ but v±[z,H]. 

The entrained steady state in terms of Eq. (2) implies 

{dSi/dzi)^ = o. (io) 

The detailed balance defined as vanishing irreversible 
probability currents between any two states of the sys- 
tem complies with condition (10). Upon (10), the equi- 
partition theorem loses its force beyond approximation 

(3) , e.g., when jump random influences is a factor, for 
the steady solution to (2) becomes non-Boltzmann, but 
[H,p st (z)] always equals zero, hence, the yardstick (1). 

The theorem IV: The canonical invariance of the 
irreversible operator of kinetic equations 

Let us show that the irreversible operator of kinetic 
equations is an invariant of canonical transformations. 
This theorem is the essence of the claimed energy dualism 
and embodies its basic rule. 

Consider a canonical transformation z — > Z of Eq. (2). 
The Poisson bracket is then to be invariant, as well as the 
probability density p(z,t). The latter turns into U(Z,t) 
according to p(z,t)d 2n z = H(Z,t)d 2n Z which results in 
p(z,t) — Yl((Z(z,t),t) since the Jacobian of any canoni- 
cal transformation 

\deb{dZ k {z,t)/dz i }% =1 \ = l. (11) 

Under the transformation, the dp/dt of (2) gives rise to 
dU(Z,t)/dt plus the addition Z' l dn(Z,t)/dZ l where 

Z' = dZ(z,t)/dt 

is the local velocity of map z — > Z. By virtue of 
(11) this velocity as a function of Z is divergence-free, 
dZ' i (z(Z,t),t)/dZi = 0, as can be shown by differentiat- 
ing both parts of (11) and using Cramer's rule. Therefore 
the addition contributes entirely to the reversible drift, 
reducing to the Poisson bracket [G, n] where 

dG(Z,t)/dX = -P', 8G(Z,t)/dP = X', (12) 

(X, P) denote the conjugated n-component coordinates 
and moments of Z, and (X', P') = Z' . That is, 

dp(z,t)/dt-[H,p] =d!l(Z,t)/dt- [H,n] 

with both sides understood as functions of either z, t or 
Z, t and where H is the transformed Hamiltonian, 

U = H + G. (13) 

So, nothing is added to the dSi/dzi part of (2) under the 
transformation and this part preserves invariance. This 
and that p is canonically invariant completes the proof of 



the theorem if S is purely irreversible, i.e. for the Hamil- 
tonian H implied dressed. The proof holds also for the 
part dSi I dzi decomposable into a pure irreversible term 
plus Poisson bracket terms. The latter are canonically 
invariant, so is therefore the irreversible term and hence 
its operator. 

It follows from the theorem that the time-symmetry 
decomposition once adopted in (2) must hold further for 
an arbitrary time dependence of the parameters of H and 
S. The irreversible operator, being invariant under ar- 
bitrary canonical transformations and initially adopted 
related to the system's motion given by the invariant 
Poisson bracket [z, H], must preserve its dependence on 
H , while the reversible operator must change as H trans- 
forms into H. This means the energy duality. The dif- 
ference G between the two Hamiltonians depends on the 
Z' given by (12) and can much exceed H in effect. Re- 
spectively the effects of energy duality can be strong. 

In general the decomposition of S is in question. Say, 
for S representing a retarded functional of p, i.e., acting 
not only in argument z but also in t, the notion of irre- 
versible operator is vague and the theorem loses its force. 
As for the energy duality principle, it is to be viewed ap- 
plicable in compliance with theorem III to an initially 
adopted H versus the operator of whole S. 



Averaging versus canonical transformations 

The emergence of irreversible behaviors is often re- 
ferred to the averaging of a conservative many-body sys- 
tem given by a microscopic Hamiltonian and random ini- 
tial conditions. This cue misleads in the question of both 
statistical and dynamical averaging. Averaging the Liou- 
vttlc's equation of the system over the statistics of initial 
conditions specifies the initial distribution function while 
its probabilistic evolution is governed by the same equa- 
tion. Also, so long as its solution is unique, integrat- 
ing this equation over the irrelevant variables reduces 
to a canonical transformation splitting the microscopic 
Hamiltonian into an averaged part independent of the 
remaining part, which leaves no place for the irreversible 
probability currents. 

While splitting of a Hamiltonian system into an av- 
eraged subsystem independent of the rest corresponds, 
ideally, to some canonical transformation, thinking the 
same way of systems where irreversible flows matter is 
incorrect. No questions arise for the averaging over the 
interactions and influences treated as small, disregard- 
ing their nonlinear effect; but the averaging causing ap- 
preciable irreversible flows is different, and it never re- 
duces to canonical transformations, for by virtue of the- 
orem IV such flows would not emerge or change. This 
is so whether the system is conservative or not. As for 
the systems under time dependent influences treated as 
a random process, its exact statistical averaging causes 
irreversible flows as set already on the probability evo- 
lution of the random process. Some way or other, the 
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irreversible flows emerge or change only due to trunca- 
tions of interactions. Just as the flows, the energy duality 
is a relative category. Regularly recurring processes lend 
weight to them, and so is the scientific cognition. 

Practically, a theory established for some conditions 
is extended further to interactions with more fields to 
model transitions to underspecified states and instabili- 
ties. The trend of irreversible flows then may change, act 
as inverse truncation if left intact, and, being uncertain, 
it is usually assumed in accordance with the paradigm 
of entrainment ideal. For example, the impact of fields 
associated with Feshbach resonance [5] on many body 
systems is commonly treated presently as the initial en- 
trainment ideal p s t — p s t(H) turning into the final ideal 
Pst = Pst(H) where Hamiltonian H includes the field. De- 
parture from this takcn-for-granted trend in transitions 
means admitting a steady non-entrained state, a relax- 
ation sort of persistent currents setting in. 

An important field where at stake is whether to stick 
to the generalized thermodynamic potential or follow the 
principles of energy duality is the impact of high fre- 
quency fields on the system's probability evolution. So 
we can view also the foregoing persistent currents in the 
moving picture where they look frozen and where the 
arising hf field is to be treated self-consistently. Let us 
dwell here on typical models of H containing hf terms 
treated as given and S of form (3) with the constraints 
on d and / as stipulated in the paragraph with Eq. (9) for 
such H = H(z, t). The hf part of H is assumed of finite 
amplitude and frequencies high enough to invalidate the 
approximation of distribution p(z,t) by quasistationary 
Boltzmann distribution (6). A characteristic measure of- 
ten in use then, and we question its extensive use, is the 
Boltzmann p(z,t) = U B [Z(z,t)], 



n B (z) 



-/3H(Z) 



e -m(Z) d 2n z 



(14) 



in terms of canonical Z where H(Z) is the quasienergy, 
the i-independent function H given by (13), or its surro- 
gate given by the effective potential i.e. the Hamiltonian 
of smoothed dynamics Z(t) obtained on hf averaging the 
equations dz/dt — [z,H(z,t)] for z — z(t). 

The conversion to terms of quasienergy or effective 
Hamiltonian is a highway in physics, it clears up the 
energy considerations of purely dynamic systems since 
the t-dependence of H makes their energy a vague no- 
tion. However, the measure (14), being the exact steady 
distribution for zero hf field, violates the canonical invari- 
ancc theorem and leaves no place for the energy duality 
for any finite hf fields, which suggests the inadequacy of 
(14) in many conditions. 

Obviously the approximation (14) is not valid on the 
stages before thermalization unless the state is prepared 
so initially. But even prepared so, the dynamic system 
Z(t) given by the Hamiltonian H(Z) is not entrained. As 
it vibrates with respect to the environment persistently 
and rapidly, in time with the hf field, there arise mean 
vortex forces having a cumulative bearing on the system's 



behaviors which is generally substantial, as elucidated 
below. 



The vortex impact of time dependent fields 

The division of generalized forces and impacts into 
types implies reasoning and transforms in terms of spe- 
cific variables. Canonical ones give division with respect 
to irreversibility, and the terms where the system Hamil- 
tonian is independent of time give a subdivision - of re- 
versible forces into potential versus gyroscopic and, as 
shown below, of their irreversible counterpart into vor- 
tex versus non- vortex forces. 

Transforming Eq. (2) with constraints (3), (9) where 
H = H{z 1 1) includes given (t dependent) fields to canon- 
ical variables Z where H (13) is t-independent, we obtain 



where the elements 



dZ k {z,t) dZ n {z,t) 



dzi 



dzj 



(15) 



(16) 



comprise the matrix D of diffusion. D is positive semi- 
definite like matrix d by virtue of (11). n, G and D 
are understood in (15) as functions of Z and t, with z 
expressed via Z and t. Obviously both z — > Z and Z — > z 
are one-to-one maps by virtue of (11). 

Since H of (15) is independent of t, the components 



dZ k 



(17) 



(where i,k = 1, ...2n) comprise the irreversible drift 
force F m corresponding to the entrainment ideal, and 
the irreversible drift force F nc of components 



Fr = (3D t 



dG 
~dZ~k 



(18) 



is due to the non-cntrainmcnt caused by the given field. 

Let us juxtapose the work done by these two forces and 
diffusion on the system. Multiplying both sides of (15) 
by H and integrating over space Z results in 

= -WD^m) + P{DikViZ< k ) + <^g*) (19) 

where V — dR/dZ is the vector orthogonal to [Z, H] and 
its norm |V"| = \[Z,U]\, and dG/dZ of (12) is rewritten 
as Z' since it is orthogonal to Z' and \Z'\ = \Z'\. 

The quadratic form D ik ViVk > 0, hence the power of 
force F el (17) is always dissipativc, tends to decrease (H). 
The mean power of diffusion forces, the third sum in the 
rhs of (19), is always positive, anti-dissipative near steady 
states in the limit of weak given fields, for only it then 
remains to balance dissipation due to F ei . 
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In contrast, the power of force F nc can be positive as 
well as negative. It is positive when the vectors V and Z' 
are parallel and negative when antiparallel. So the sign of 
work of F no over a closed path in phase space Z depends 
on direction of path-tracing, which is intrinsic of a mean 
vortex force. As evident from (19), the mean power of 
vortex force prevails over that of dissipative F el as the 
speed V of system decreases compared to the speed Z' 
of non-entrainment. 

Not only the steady distribution can then strongly dif- 
fer from Boltzmann (14), but also the stability threshold 
of the system can appreciably retreat toward as inside as 
outside the domain given by (14), i.e., where the H(Z) is 
bounded below. Also the rates of relaxation can be much 
higher or lower than given by eigenvalues of [3D. 

Particularly this concerns the nonlinear impact of high 
frequency fields, even relatively weak, as is the case 
of long lasting rotation/vibration systems near reso- 
nances, including parametric and combinational. This 
was brought up in various fields by the present author, 
e.g. ,[6-12]. Note that we referred to as vortex there only 
the vortex force field in the coordinate space of system, 
rather than in all phase space, as in the present work. 



Multi-bath extension 

The energy duality formulated above is with respect 
to the guide of entrainment limit. It does not cover the 
realm of environments not carried along between them- 
selves on the average. However, the energy duality prin- 
ciple then displays through treating the realm as a super- 
position of environments ("baths") each compliant with 
detailed balance and moving with respect to each other. 
Applied to diffusional approximation, this is a number 
of baths each compliant with constraints (9) on d and / 
in the picture of canonical variables where the system is 
carried along it. On account of the canonical invariance 
theorem we then come out with the irreversible drift force 
/ related to the diffusion matrix d in Eq. (3) so 



(20) 



Here means summing over the superscript r labeling 
the quantities related to each of the baths, /3 r is the noise- 
energy scale of bath r, d r its contribution to the diffusion 
rates of the system, and H + G r the Hamiltonian trans- 
formed to the picture of canonical variables z r at rest 
with bath r. With G r given in this picture as function of 
z r and t, we have 

<9G r (z r (z,i),t) _ dzl(z,t)dzl(z,t) 



dz k 



dz k 



dt 



(21) 



where z r = (— p r ,x r ), x r are the coordinates and p 1 the 
moments of z T = (x r ,p r ). 

It follows from (20) that the components 



dH{z,t) 
dz k 



comprise the irreversible drift force / cl exerted on the 
system in neglect of its motion relative to the baths, while 
the force f ne of components 



fi 



(23) 



with dG v /dzk given by (21) is due to the relative motion 
of the system. The / nc is vortex in space z and / el is non- 
vortex, similar to the hf forces (18) and (17) in space 
Z. We now arrive in space Z with the energy balance 
equation 



<*(H) 
dt 



J2P r (Dl k V t (Z' k 



Z'k)) 



(22) 



where D\ k is of form (16) with dkn replaced by d kn and 
Z n stands for dG'/dZ = JdG r /dz with J the inverse of 
matrix dz(Z, t)/dZ. 

As with (19), the first sum on the right of (24) is always 
dissipative, tends to decrease the mean energy (H) of the 
system and the second sum represents the power of the 
mean vortex force whose work over a closed path in phase 
space Z depends on direction of path-tracing. 

Remarkably, the relation (20) between the irreversible 
drift force and diffusion for dG r /dz = for all baths 
looks like the constraint (9) with f3 = ^/3 r d _1 d r . But, 
this [3 is a matrix function of z rather than a constant 
scalar. It is easy to show that the relation (9) with such 
(3 represents a general constraint for the systems in the 
entrainment limit of diffusional approximation. 

It follows that the equipartition theorem in such con- 
ditions of entrainment generally does not take place and 
that the steady distribution p st can be reduced to a super- 
position of a number (< 2n) of Boltzmann distributions. 
The latter form of solution breaks beyond the entrain- 
ment limit as well, which can be used in determining the 
limit. 

To illustrate features of the energy duality and asso- 
ciated vortex physics arising in non-entrainment condi- 
tions, a simple example relevant to applications is con- 
sidered below. 



The vortex confinement of particles in fields of 
rotational symmetry 

Here we apply our approach to long time confinement 
of particles with fields that are constant in some rotating 
frame. For simplicity the particle trap formed in axial 
direction by the field is assumed harmonic and indepen- 
dent of transversal particle motion. The transversal force 
field is allowed of arbitrary azimuthal asymmetry. So are 
Penning and many other traps for charged particles and 
neutral atoms, including a general quadrupole type of 
traps suggested by the author in 90-ties, see references 
in [12] devoted to elemental vortex statistical mechanics 
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of relevance. While the theory of generalized thermody- 
namic potential gained ground in such traps [13], it seems 
to be severely restricted in view of non-entrained back- 
ground gas, noisy field sources and retarded reaction of 
trap field system (due to its finite conductance) to the 
particle motion. 

Let us rely on [12], bringing in a wider insight and 
dwelling on the probability of particle motion governed 
by Eqs. (2), (3) with constraints of form (4) for a number 
of baths each rotating about the symmetry axis of trap 
at some frequency ft r and of constant, isotropic diffusion 
rates d r . In the rotating- frame coordinates where the 
trap field is independent of t, the kinetic equation (15) 
for the probability distribution of particle states reduces 
to 

|n(x, p, t) = [h, n] + A (p DVi ft* + D-^j n 

'(25) 

where 

D = Y,<? and /3 = l^/3^. 



D 

V = dU/dP, —(3DV is the friction force F oi (17), F nc 
(18) takes the form of net vortex force F vort in space X. 
For a particle of Euclidean coordinates X = (X\,X 2 ,X 3 ) 

F vmt = (kX 2 , -kX u 0), k = ^(fi - ft r )/3 r d r , (26) 

ft is the rotation frequency and X 3 along the rotation axis 
of trap. In the trap field modeled as given, the Hamilto- 
nian of particle motion is of form 



H = ^[P 2 + k 3 X 2 3 + (ko-2q)X 2 1 +(ko + 2q)Xl] 



+ (<?-ft)M e + H 



ni 



(27) 



where P = (P\,P 2 ,P 3 ) are canonical moments conju- 
gated to X. X,P are related to the canonical x,pof rest 
frame so X = T(t)x, P = T(t)p, 



T(t) 




sin Clt 
cos fit 
1 



G = H — H = —CIMq where H is the rest- frame Hamil- 
tonian and Mq the canonical angular momentum 



Me = xip 2 - x 2 pi = X\P 2 - X 2 Pi, 



(28) 



ko = g 2 — k 3 /2, k 3 > 0, y/k~3 is the eigenfrequency of ax- 
ial oscillator mode, q the gradient of rotating transversal 
quadrupole potential and H nl the potential of transver- 
sal multipoles higher than quadrupole. H nl = for the 
geometries of Penning and rotating quadrupole traps. 

The confinement in line with the theory of general- 
ized thermodynamic potential corresponds to the ideal 
of k = and would mean that the particle should relax 
to the entrained state of rotating-frame Boltzmann dis- 
tribution (14) in the parameter domain where H(Z) (27) 



is bounded below, and an instability should arise as the 
boundedness breaks down; for H nl = 0, it breaks down 
where 



k < (ft - gf 



(29) 



The critical point ko = (ft — g) 2 is the apex of paramet- 
ric resonance of transversal motion caused by the rotat- 
ing quadrupole. However, the trends of long time con- 
finement in considered traps strongly diverge from these 
predictions. For example, the traps of g — would be 
impossible in principle, for any q, if condition (29) would 
indeed imply instability, but it may not at all in view of 
the non-entrainment given by n 7^ 0. 

An important point to account for k 7^ is that the 
work of vortex force over a closed path enclosing in a 
transversal plane (across axis X 3 ) an area s is equal 
to § F vmt dX — ±2ks, velocity independent, and the 
sign depends on the direction of motion, while the work 
§ (3DVdX decreases with V and is of one sign. For a 
circular transversal motion of a frequency u) 



§F vmt dX 1 _ _ Q* 
§ -f3DVdX ± ~ ~ 



where X± denotes a vector (X\,X 2 )- A sizable propor- 
tion of CI* is due to the environments roughly at rest with 
the trap electrodes, of £! r <C O, so CI* and CI are within 
one order. The soft mode of transversal dynamics given 
by H is slow, of frequency scale v <C CI. So ui <C CI* for 
such modes. The particle stability, relaxation and steady 
state is critical to the soft mode, hence the dominance of 
vortex force impact on all that. 

Importantly the dominance holds on passing the crit- 
ical point into a certain parameter domain of H(Z) un- 
bounded from below, where the soft mode still represents 
a transversal orbiting, rather than runaway motion. Such 
a domain is characteristic of dynamic systems with gyro- 
scopic terms, as is the case of particle traps with rotat- 
ing fields. For H nl = 0, the Hamiltonian of transversal 
dynamics presented in normal mode presentation in the 
domain takes the form [12] 

Hl = \v + {P% + X 2 + ) ± \v_(P 2 _ + X 2 ) (31) 

where '±' is '+' for k > (Cl-g) 2 and '-' for k < (Cl-g) 2 , 
X + ,P + arc canonical variables of normal '+' mode and 
X-, P- of ' — ' mode, v± are positive roots of 



4 = fc + (ft - gf ± 2V(ft-5) 2 fco+<Z 2 , (32) 

and the domain of no runaway is where 

v 2 _ > and k 3 > 0. (33) 

For the apex falling in domain (33), the '— ' mode is soft 
near the apex and its energy is negative below it. 

Along with the dominance, another important feature 
of vortex force displays vigorously and not alike the fric- 
tion force — (3DV . The latter pumps energy out of any 
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motion, so its pumping out of positive-energy modes 
causes their damping, and pumping out of negative- 
energy modes causes oscillation build-up, a negative fric- 
tion. So, on accounting only for the fiction force, the con- 
dition (29) would imply instability of increment ~ (3D. 

In contrast, the stabilizing/destabilizing impact of the 
vortex force does not change sign on passing the critical 
point but changes it on passing the point of gyroscopy 
compensation g = fl by field rotation, unless both points 
coincide. Indeed, the energy pumped out/in the soft 
mode by the vortex force depends on the ellipticity of 
transversal-mode orbiting and whether in the rotating 
frame the orbiting is in the direction of field rotation or 
opposite to it. For g ^ ft, the sense of orbiting changes 
twice, on passing both points, but on passing the crit- 
ical point there changes the sign of energy transfer by 
the vortex force, hence, its stabilizing/destabilizing im- 
pact does not change sign on transition into the domain 
of H(Z) not bound below for the same reason why there 
changes the sign of friction force impact. 

For the models considered in this section, the vortex 
impact vanishes for the transversal modes of linear polar- 
ization, which is when g — Q, and is maximal for circular 
orbiting. For H nl = 0, calculating the time averages of 
energy flows due to the vortex and friction forces in the 
orbiting motion corresponding to the '+' and '— ' modes 
with exact account of their polarization forms, we obtain 
for their ratio, respectively 

/ l F^dx^ \ n* (g - m (U) 
\§-(3DVdx ± J ± n 4-feo ' 1 ' 

The ratio O* /fl > 0, so the sign of ratio next to fi*/f2 in 
(34) determines the sign of stabilization/destabilization 
impact of the vortex force. There the dominator v\ — fco 
changes sign only on passing the apex for the '— ' mode, 
as evident from (32). Thus the expressions (34) and (33) 
show all variety of confinement trends in question. 

For <? = 0, for example, fco < and the domain (33) is 
completely on the side of apex where (29) holds and the 
soft mode is of negative energy. However, the vortex force 
exerts there stabilizing action on the soft mode, and for 
0* = O it prevails over destabilizing effect of force —(3DV 
in all domain (33). As for the '+' mode, the vortex 
impact on it is destabilizing but weaker than frictional 
which is stabilizing. So, a stable 3D confinement takes 
place in all domain (33)! 

At the point of gyro-compensation g = fl, the vor- 
tex factor (34) vanishes and the trends of confinement 
in the trap fields of H nl = comply with the theory of 
generalized thermodynamic potential. However, this spe- 
cific point is aside from conditions most favoring particle 
cooling and confinement in traps, as well as heating and 
selection of particles, and there the vortex force rules the 
trends. Depending on the ratio g/tt and its sign, the 
vortex impact appears to be either stabilizing or destabi- 
lizing, resulting in amazing trends of stability, relaxation 
rates and steady states. 



Highly important in all that is the scale f2*/f2, whether 
it is small enough to stick to the entrainment-ideal sce- 
nario or large, comparable with 1, to have an essen- 
tially vortex picture of confinement. Characteristically 
the scale and hence the factor of prevalence of the 

vortex impact over frictional is the same for any positive 
values of baths' parameters (3 and D, including the limit 
of (3D -> 0+. 

For H nl = 0, the kinetic equation (25) admits exact 
analytical solution H(Z, t). It represents, for the system 
evolution from a given initial state in a point of phase 
space Z, a multivariate Gaussian distribution of system 
states for any t > 0. So all probability characteristics of 
transient and steady states can be traced exactly for the 
entire range from = to 1. 

For H nl the transversal Hamiltonian as nonlinear is 
generally inseparable into normal modes, but in any close 
parameter vicinity of inseparable Hamiltonian there exist 
separable ones. So for fairly long times the motion in 
the domain of no runaways can be treated as separable 
and of soft mode on approaching apexes of parametric 
or combination resonances. Therefore the vortex trends 
resembling those shown above are to be expected. 

Additional remarks. For more complex objects than 
a 3D-particle in a field of rotational symmetry, represent- 
ing nonlinearly coupled systems having many degrees of 
freedom of high and low frequencies, the vortex impact 
appears to be no less significant, as is, e.g., in condi- 
tions of hf subsystems excited near main or parametric 
or combinational resonances. Its general features in the 
lowest order of nonlinearities were elucidated in [6-12] 
where various methods of cooling and control of inter- 
acting particles, waves, domain walls etc were suggested 
on this basis. 

Imposing on a system in conditions of entrainment a 
regularly varying field modeled via Hamiltonian interac- 
tions always translates into a vortex impact exerted on 
the system in the picture where the field looks frozen. 
This refers to any control of irreversible effects this way. 
So works the principle of energy duality. Now, with this 
principle, there arises a consistent and adequate energy 
approach to the above-cited and many other phenomena 
in open systems. 



Conclusion 

We have shown that the notion of energy function 
developed in analytical mechanics and thermodynamics 
acquires rigor and significance not only in the limit of 
detailed balance but also in the statistical approach to 
the vast world where it rules relaxation to non-entrained 
states, and that this comes with breaking the energy 
function into two having different transformation rules. 

The two functions merge in the entrainment ideal and 
can strongly differ beyond it. As demonstrated, both 
are important for the energy approach to the long time 
behavior and stability of systems. Otherwise, huge in- 
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consistencies arise, as is the case of particle confinement 
and other systems mentioned in this work. 

Behind the energy duality there stands the persistent 
irreversible drift related to the difference between the 
two energy functions of system and representing a vor- 
tex force field in its phase space. It is essential that the 
vortex force effect differ radically from that of familiar 
friction of systems relaxing to entrained states. We see 
its distinctive features as in the introduced classification 
of generalized vortex forces and energy flows as in the 
physics the approach gives to various phenomena. 

The principle of energy duality does not rely on quan- 
tum mechanics and one cannot but infer from the corre- 



spondence principle that a consistent concept of energy 
quanta is to be either limited by the entrainment ideal 
or the energy quanta break into two respective sorts be- 
yond the ideal. This disputes common reasonings of non- 
entrained states in terms of energy quanta. 

While the generalized thermodynamic potential theory 
fits only the entrainment limit, it may look all-sufficient 
beyond it with fitting via extending the set of variables 
and parameters. However, the theory then loses its pre- 
dictive strength and causes misleading associations like 
the ones revealed in this work. The energy duality is 
free of these drawbacks and presents a universal energy 
concept incorporating the vortex physics. 
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